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Abstract 

A  computational  model  is  developed  to  describe  the  fluid  flow  in  the  small  channel  network  systems  that  are  typical  of  planar  fuel  cells.  The 
model  solves  mass  and  momentum  equations  to  predict  pressure  drop  and  flow  uniformity  within  individual  channels.  By  formulating  the 
problem  in  dimensionless  variables,  the  solutions  are  generalized  in  terms  of  two  nondimensional  groups.  The  results  are  presented  in 
graphical  form,  facilitating  their  application  without  the  need  to  solve  the  model  for  any  particular  application.  ©  2002  Elsevier  Science  B.V. 
All  rights  reserved. 
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1.  Introduction 

Fuel  cells  are  emerging  as  a  practical  technology  for  high- 
efficiency  electric-energy  generation.  A  planar  architecture 
provides  an  effective  approach,  achieving  high  utilization  of 
the  electrochemically  active  area.  Small  channels,  which  are 
often  formed  in  the  interconnect  material,  carry  the  fuel  and 
air  gas  flow.  A  designer  must  decide  how  to  configure  a 
channel  network,  such  that  the  appropriate  flow  rates  and 
flow  uniformity  can  be  maintained  across  the  surfaces  of  the 
membrane-electrode  assembly  (MEA).  An  important  con¬ 
sideration  is  the  net  pressure  drop,  which  should  be  as  low 
as  possible  to  reduce  parasitic  power  needed  to  drive  pumps 
or  compressors. 

There  is  a  great  deal  of  current  research  on  fuel  cells 
generally,  especially  considering  new  material  systems  and 
electrochemical  performance  of  the  MEA  assembly.  How¬ 
ever,  relatively  little  of  the  literature  concentrates  on  the 
system-level  layout  of  the  layers  and  stacks.  Iwata  et  al.  [1] 
provides  one  of  the  best  overall  analyses  of  a  planar  solid- 
oxide  fuel-cell  (SOFC)  system.  Boersma  and  Sammes  [2] 
describe  a  model  to  predict  flow  in  fuel-cell  channel  net¬ 
works.  Their  model  is  based  on  networks  of  hydraulic 
resistances.  However,  the  results  of  the  model  are  not 
generalized  or  reported.  Therefore,  to  apply  the  approach, 
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one  needs  either  to  write  a  program  that  implements  the 
model  or  obtain  the  authors’  software. 

The  objective  of  this  paper  is  to  assist  the  designer  in 
making  quantitative  decisions  for  the  various  channel 
dimensions.  The  analysis  is  based  on  a  computational  model 
that  solves  mass-  and  momentum-conservation  equations 
within  the  channel  network.  Formulating  the  model  in  a 
dimensionless  framework  enables  a  generalization  of  the 
results  in  terms  of  two  dimensionless  parameter  groups.  The 
nondimensional  results  are  presented  in  graphical  form, 
establishing  practical  guidelines  for  channel  layout  and 
sizing. 

While  the  principal  intent  of  this  work  is  to  assist  the 
design  of  fuel-cell  systems,  it  should  be  noted  that  other 
applications  can  benefit  from  the  results.  For  example,  plate- 
and-frame  heat  exchangers  often  have  a  header-channel 
configuration  that  is  similar  to  a  planar  fuel-cell  layer.  In 
that  case,  the  hot  and  cold  fluid  channels  take  the  place  of  the 
fuel  and  air  channels  in  the  fuel  cell. 

1.1.  Planar  SOFC  architecture 

Fig.  1  illustrates  the  principal  features  of  a  planar  SOFC 
architecture.  The  membrane-electrode  assembly  (MEA) 
layer  is  composed  of  an  electrolyte  that  is  sandwiched 
between  electrodes.  The  MEA  layer  is  typically  on  the  order 
of  1  mm  thick,  the  details  varying  from  one  design  to 
another.  The  electrodes  (anode  and  cathode)  are  porous 
cermet  materials  that  permit  gases  to  be  transported  between 
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Fuel  channel 


Rib 


1.2.  Assumptions 


Fig.  1.  Illustration  of  a  segment  of  two  fuel-cell  layers,  including  the 
electrode-electrolyte  assembly  (MEA),  the  interconnect,  and  the  flow 
channels  for  fuel  and  air.  A  counter-flow  situation  is  illustrated  here,  but 
co-flow  and  cross-flow  configurations  are  also  common. 


the  flow  channels  and  the  electrolyte  surface.  The  electrodes 
are  loaded  with  a  metal  that  serves  two  purposes.  One  is  to 
catalyze  electrochemical  reaction  at  the  electrode-electro¬ 
lyte  interfaces,  and  possibly  to  promote  catalytic  reaction  of 
the  fuel  (e.g.  reformation  of  hydrocarbon  and  water  vapor  to 
carbon  monoxide  and  hydrogen).  The  other  purpose  of  the 
electrode  loading  is  to  carry  electric  current  from  the 
electrolyte  to  the  metal  interconnect  structure.  As  shown 
in  Fig.  1,  the  flow  channels  are  formed  in  the  interconnect 
metal. 

The  ribs,  which  separate  the  flow  channels,  make  direct 
contact  with  the  electrodes.  In  designing  the  layer  architec¬ 
ture,  there  is  a  tradeoff  that  must  be  considered  between  the 
rib  and  channel  sizes.  Thicker  ribs  may  reduce  the  interface 
resistance  to  current  flow  by  increasing  the  electrode-inter¬ 
connect  contact  area  and  by  decreasing  the  current  path 
through  the  relatively  high  resistance  electrode  material. 
From  a  fluid- flow  point  of  view,  wide  channels  and  narrow 
ribs  facilitate  more  uniform  distribution  of  reactive  gases 
across  the  full  area  of  the  electrolyte  and  thus,  promote 
electrochemical  performance.  The  results  of  this  paper  assist 
the  designer  in  understanding  how  channel  sizing  affects 
flow  uniformity  and  pressure  drops  for  the  system. 

A  single  MEA  layer  typically  produces  an  areal  power 
density  that  is  on  the  order  of  1  W/cm2.  Open  circuit 
voltages  are  on  the  order  of  1  V,  with  maximum  power 
output  at  roughly  1/2  V.  Maximum  current  densities  are 
typically  on  the  order  of  5  A/cm2.  The  performance  of 
any  particular  system  varies  depending  on  a  great  many 
design  and  operational  considerations.  Because  the  layer 
voltage  is  less  than  1  V,  it  is  usually  necessary  to  connect  the 
layers  in  series  to  increase  the  voltage  of  a  stack.  As 
illustrated  in  Fig.  1,  layers  can  be  stacked  one  on  top  of 
another  much  as  stacking  batteries  in  a  flashlight.  To  deliver 
a  stack  voltage  of  say  24  V,  over  50  layers  are  needed.  Thus, 
in  addition  to  understanding  flow  distribution  within  a  layer, 
one  must  also  design  to  produce  uniform  flow  to  each  layer 
in  the  stack. 


The  analysis  here  presumes  an  incompressible  flow,  i.e. 
the  mass  density  p  is  constant.  While  the  fluid  is  typically  a 
gas  (e.g.  air  or  vaporized  fuel),  the  flow  velocities  are  quite 
low  and  the  temperature  is  nearly  uniform.  Moreover,  all 
other  fluid  properties  (e.g.  viscosity)  are  taken  as  constant. 
While  incorporating  variable  properties  presents  no  princi¬ 
ple  problems  in  the  computational  analysis,  using  constant 
properties  facilitates  a  generalization  of  the  results  via 
characteristic  nondimensional  groups. 

The  following  sections  of  this  paper  develop  the  govern¬ 
ing  equations  and  discuss  their  computational  solution  from 
the  point  of  view  of  flow  across  a  single  fuel-cell  layer. 
Flowever,  there  is  also  a  need  to  consider  the  flow  through  a 
stack,  which  is  composed  of  many  layers.  Fortunately,  the 
stack  analysis  follows  directly  from  the  layer  analysis.  For 
the  stack  analysis,  the  layer  takes  the  role  of  the  channels  in 
the  layer  and  the  manifolds  take  the  role  of  the  headers. 

The  analysis  here  considers  flow  in  impervious  channels. 
It  does  not  account  for  the  mass  flow  through  the  electrolyte, 
and  thus,  the  mass  flow  between  the  oxidizer  and  fuel 
channels.  Certainly,  this  is  a  significant  assumption  for  a 
fuel-cell  application,  the  operation  of  which  depends  on  the 
ion  flux  between  the  fuel  and  air  channels.  Nevertheless,  the 
analysis  represents  a  conservative  limiting  case  for  sizing 
the  air  channels  (which  lose  mass  in  an  SOFC  cell).  Because 
the  air-flow  rates  are  so  large  compared  to  the  fuel  flow,  it 
is  the  air  channels  that  typically  demand  the  most  design 
attention  to  achieve  uniform  flow.  The  intent  of  this  analysis 
is  to  provide  channel-sizing  guidance  in  a  generalized  form. 
Incorporation  of  electrochemistry  and  the  consequent  mass 
exchange  frustrates  the  ability  to  generalize  the  results,  and 
thus,  provide  easy-to-use  design  guidance. 


2.  Governing  equations 

The  flow  distribution  across  a  layer  (e.g.  Fig.  2)  is  based 
on  solving  mass-continuity  and  momentum  equations  along 
the  length  of  the  feed  and  exhaust  headers.  Fluid  flows 
between  the  feed  and  exhaust  headers  via  the  channels  that 
connect  them.  Thus,  the  conservation  equations  in  the  head¬ 
ers  must  accommodate  the  distributed  mass  loss  from  the 
feed  header  and  the  mass  increase  in  the  exhaust  header. 

It  is  clear  from  Fig.  2  that  there  are  a  certain  number  of 
channels  in  the  layer  and  from  Fig.  3  that  there  are  a  certain 
number  of  layers  in  the  stack.  With  this,  the  system  of 
equations  are  sensibly  formulated  and  solved  as  a  system  of 
discrete  algebraic  equations  that  align  with  the  channels  or 
the  layers.  However,  we  find  it  beneficial  to  first  derive  the 
equations  as  systems  of  continuous  differential  equations.  In 
a  differential-equation  form,  it  is  easier  to  understand  the 
mathematical  structure  and  the  requirements  for  boundary 
conditions.  The  solution,  however,  is  accomplished  in  discrete 
form. 
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Fig.  2.  Illustration  of  the  flow  distribution  in  a  fuel-cell  layer,  showing  the  manifolds  that  carry  flow  between  layers  and  the  headers  that  distribute  flow  to  the 
individual  channels  across  the  layer. 


The  analysis  in  the  following  sections  concentrates  on 
the  layer.  An  analogy  is  then  drawn  between  the  layer  and 
the  stack,  showing  how  the  layer  analysis  is  generalized 
to  the  stack. 


can  be  represented  in  terms  of  a  nondimensional  friction 
factor  /  as 


(2) 


2.1.  Hagen-Poiseuille  channel  flow 

For  steady  Hagen-Poiseuille  flow  in  a  channel,  a  net 
pressure  variation  is  needed  to  balance  the  force  caused 
by  the  shear  stress  at  the  walls 

( Pi  Pc  Me  =  rw/JcLc  (1) 

where  pf  and  pe  represent  pressures  in  the  feed  and  exhaust 
headers,  i.e.  the  inlet  and  outlet  of  a  channel.  Other  variables 
are  Ac,  the  channel  cross-sectional  area;  Pc,  the  channel 
perimeter;  and  Lc,  the  length  of  the  channel.  The  shear  stress 


Fig.  3.  Illustration  of  a  full  stack,  with  many  individual  layers.  The  fuel 
and  air  manifolds  must  feed  each  layer,  and  collect  the  exhaust  from  each 
layer. 


where  u  is  the  mean  velocity  in  the  channel.  For  steady 
laminar  flow  in  a  channel,  the  product  of  the  friction  factor 
and  the  Reynolds  number  (based  on  mean  velocity  and  a 
hydraulic  diameter)  is  a  constant. 

Ref  =  constant  (3) 


where,  for  a  circular  channel  Ref  =16.  The  Reynolds 
number  and  the  hydraulic  diameter  are  written  as 


Re  = 


puDc 

P 


(4) 


For  rectangular  channels,  the  value  of  Ref  depends  on  the 
channel  aspect  ratio  a  =  L/W  as 


Ref 


13.74+  10.38  exp 


(5) 


In  a  square  channel  Ref  «  14,  while  for  a  very  high  aspect- 
ratio  channel  Ref  ss  24  [3], 

Yuan,  et  al.  [4]  have  simulated  the  flow  in  variously 
shaped  fuel-cell  channels.  They  report  both  Ref  values  as 
well  as  Nusselt  numbers  for  heat-transfer  analysis.  In  addi¬ 
tion  to  studying  channel  shapes,  they  have  also  characterized 
the  effects  of  mass  transfer  through  the  channel  walls  as 
happens  in  a  fuel-cell  channel.  The  results  show  that  the 
mass  transfer  through  the  channel  walls  has  a  relatively 
small  influence  on  the  value  of  Ref. 

With  a  certain  amount  of  algebra,  an  expression  can  be 
developed  to  relate  the  mass-flow  rate  in  the  channel  to  the 
net  pressure  difference  from  the  channel  inlet  to  the  outlet 


mc 


(  1  d2Acp\ 
\Ref  2 Lc  p) 


Of  ~  Pe) 


Kc(j>f  -pe) 


(6) 
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where  /<  is  the  dynamic  viscosity  of  the  fluid.  The  constant  Kc, 


1  P2cAcp 
Re/  2 Lc  / i 


(7) 


characterizes  the  proportionality  between  the  mass-flow  rate 
in  a  channel  and  the  net  pressure  difference  from  inlet  to 
outlet.  In  the  fuel-cell  layer,  pf(z)  and  pe(z)  are  both  func¬ 
tions  of  position  in  the  header.  Therefore,  depending  on  the 
pressure  distributions  in  the  headers,  the  mass  flow  in  each 
channel  may  be  different. 

In  the  foregoing  analysis  the  flow  conductance  Kc  is 
attributed  to  viscous  effects  in  the  channels.  However,  as 
discussed  by  Boersma  and  Sammes  [2],  other  factors  can 
contribute  to  flow  resistance.  These  include  expansion  and 
contraction  as  the  flow  enters  or  leaves  a  channel  or  the 
resistance  associated  with  the  gas  changing  direction  as  it 
leaves  the  header  and  enters  a  channel.  Imperfections  or 
roughness  in  the  channel  walls  can  also  increase  flow 
resistance.  All  these  effects  can  be  incorporated  into  the 
present  analysis  by  suitably  reducing  the  value  of  Kc. 

As  noted  earlier,  the  equations  are  first  derived  as  con¬ 
tinuous  differential  equations.  Thus,  it  is  necessary  to  repre¬ 
sent  the  mass-flow  rate  between  the  headers  as  a  continuous 
variable.  The  channel  mass-flow  rate  per  unit  length  of 
header  (kg/m-s)  may  be  written  as 


m 


I 

C 


K-ciPi-Pe) 

K 


K  (Pf-Pe) 


(8) 


where  A'  is  the  channel-opening  cross-sectional  area  per  unit 
length  of  the  header.  In  a  limiting  case,  one  can  imagine  an 
infinite  number  of  channels,  such  that  the  flow  between  the 
feed  and  exhaust  headers  varies  continuously  along  the 
length  of  the  headers. 


2.2.  Continuity  equation  in  headers 


Derivation  of  the  continuity  equation  is  formulated  for  a 
differential  control  volume  within  a  header  (Fig.  4).  Beginning 
with  the  Reynolds  Transport  Theorem  (relating  a  system  to  a 
control  volume)  and  the  general  mass-conservation  law  for  a 
system,  the  continuity  equation  is  given  as 


(9) 


where  m"  is  the  mass  flux  that  leaves  (or  enters)  the  headers 
through  the  connecting  channels.  Note  that  the  mass  flow 
through  the  channels  (i.e.  between  the  headers)  alters  the 
system  mass  within  the  header  control  volume.  For  the 
constant  cross-section  headers,  dV  =  Ahdz,  where  Ah  is 
the  cross-sectional  area  of  the  header  (for  a  rectangular 
channel  Ah  =  H^WO- 

The  area  through  which  the  channel  flow  leaves  (normal 
to  the  header)  is  dA  =  A'  dz.  As  illustrated  for  the  rectangular 
channels  in  Fig.  4,  A'  =  Hc.  Shrinking  the  control  volume  to 
a  differential  length  dz,  the  integrands  can  be  considered  as 
constants.  Dividing  by  6V  =  Ah<5z  and  restricting  attention 
to  steady-state,  one-dimensional  flow  provides  the  desired 
continuity  equation  in  differential-equation  form 


d  (pu) 
dz 


A' 


(10) 


Assuming  incompressible  flow  (i.e.  p  is  constant), 
dz  pAh  pAh 


(11) 


Substituting  Eq.  (8),  the  continuity  equations  take  the  form 


dwh  =  ±K(Pt~Pe) 
dz  pAh 


(12) 


where  the  subscript  h  =  f  refers  to  the  feed  header  (right 
hand  side  takes  the  minus  sign)  and  h  =  e  refers  to  the 
exhaust  header  (right  hand  side  is  positive). 


2.3.  Momentum  equation  in  headers 


In  general,  the  momentum  balance  is  stated  as 


where  the  extensive  variable  is  the  momentum  vector 
P  =  mV  and  F  is  a  force  exerted  on  the  system.  The 
intensive  variable  associated  with  the  momentum  is  the 
velocity  vector  V.  Thus,  for  one-dimensional  axial  flow  in 
the  headers,  the  Reynolds  Transport  Theorem  provides 


d{puu)\ 
&  ) 


(14) 


There  are  two  forces  that  need  to  be  considered.  One  is  due 
to  the  pressure  acting  on  the  cross-sectional  flow  area  and  the 
other  comes  from  viscous  drag  at  the  channel  walls.  Restrict¬ 
ing  attention  to  steady  flow, 


f  d(puu) 
./cv  dz 


■  ndA  — 


(15) 


The  first  integral  on  the  right-hand  side  represents  the  force 
due  to  normal  stress  on  the  control-volume  faces  in  the  flow 
direction.  In  this  case,  the  normal  stress  is  due  entirely  to  the 
pressure,  xz  =  —p.  The  second  integral  represents  the  shear 
drag  at  the  channel  walls.  The  Gauss  Divergence  Theorem 
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can  be  used  to  convert  the  first  integral  on  the  right-hand  side 
to  a  volume  integral.  The  second  integral  is  rewritten  in 
terms  of  the  differential  length  dz. 


lew 


d(puu) 

~lk~ 


dV  = 


rzdV  - 


TwFhdz 


lew 


(16) 


where  //  is  the  perimeter  of  the  header  channel  (for  the 
rectangular  channel  in  Fig.  4,  P h  =  2//h  +  2Vkh).  Shrinking 
the  control  volume  to  vanishingly  small  <5z,  yields  the 
desired  differential  equation  as 

d(puu)  dp  Ph 

A  —  A  A  Tw  (1  ') 

dz  dz  Ah 


For  an  incompressible  flow,  the  equation  can  be  rewritten  as 


2  pu 


du 

dz 


dp 

dz 


=  0 


(18) 


As  in  the  channel  flow,  the  wall  shear  stress  tw  is  determined 
in  terms  of  a  friction  factor/,  which  depends  on  the  header 
dimensions  and  fluid  properties. 

By  substituting  the  expression  for  shear  stress,  the 
momentum  equation  can  be  rewritten  as 


2piih 


dui, 

dz 


Ph 

2AhDh 


Mh 


0 


(19) 


Assuming  constant  properties  is  seen  that  the  shear-stress 
term  is  the  product  of  a  constant  factor  and  the  local  velocity. 
For  the  sake  of  simplifying  the  nomenclature,  the  constant  in 
the  shear-stress  term  may  be  defined  as 


Fh 


Ph 

2AhDh 


/(Re/ 


(20) 


As  in  the  continuity  equation,  a  subscript  h  has  been 
introduced  to  represent  either  the  feed  or  the  exhaust  header. 
In  the  subsequent  analysis,  it  is  presumed  that  the  feed  and 
exhaust  headers  have  the  same  dimensions,  hence  the  same 

Fh. 


Fig.  5.  Projection  view  of  the  header  and  channel  layout,  showing  the  unit 
cells  of  the  finite-volume  mesh  network.  In  this  illustration  the  headers  are 
divided  into  seven  unit  cells,  with  eight  mesh  points.  The  boundary 
conditions  are  shown  at  the  ends  of  the  headers. 


2.5.  Boundary  conditions 

Since  all  the  equations  are  first-order  ordinary  differential 
equations,  each  demands  one  boundary  condition.  Thus,  in 
total,  four  boundary  conditions  are  required.  As  indicated  in 
Fig.  5,  the  velocity  must  be  zero  at  the  closed  ends  of  the 
headers 

Uf(Lh)  =  0,  ne(0)  =  0  (25) 

The  pressures  in  the  feed  and  exhaust  manifolds  can  also 
provide  suitable  boundary  conditions. 

pf{0)=pQ,  /?e(Lh)=0  (26) 

The  exhaust-manifold  pressure  may  be  taken  as  a  reference 
pressure  of  zero.  However,  for  the  application  here,  it  is 
appropriate  to  specify  the  net  mass-flow  rate  M\  for  the  layer 
instead  of  the  pressure  in  the  feed  manifold.  Therefore,  the 
feed  pressure  becomes  an  implicit  boundary  condition  that 
must  be  determined  during  the  course  of  the  solution. 


2.4.  Equation  summary 


3.  Nondimensionalization 


In  total,  a  system  of  four,  coupled,  nonlinear,  first-order, 
ordinary  differential  equations  describe  the  header  flow 


dMf  =  KiPt-Pe) 
dz  pAh 


dMe 

dz 


IT 

i 

CD 

Mh 

d«f  dpf  j 

(  Ph 

dz 

dz  ' 

V  2 AhDh 

di/e  dpe 

(  Ph 

e  dz 

dz 

\  2AhDh 

/(Re/  J  u f 
pRef]ue 


(21) 

(22) 

(23) 

(24) 


As  illustrated  in  Fig.  5  for  a  layer,  the  independent-variable 
domain  spans  0  <  z  <  /,. 


With  the  objective  of  generalizing  the  results,  the  equa¬ 
tions  can  be  recast  in  terms  of  nondimensional  variables.  The 
reference  length  scale  is  the  header-channel  length  Lh  and 
the  reference  velocity  is  derived  from  the  layer  mass-flow 
rate  M\  as 


(27) 


The  reference  velocity  is  the  mean  velocity  that  enters  the 
feed  header  and  exits  the  exhaust  header.  The  reference 
pressure  is  a  dynamic  pressure  based  on  the  reference 
velocity 


,  M, 

Pk{  ~pU'~^ 


(28) 
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The  nondimensional  variables  are  thus  defined  as 

2  -  U  Mh  ,  p  p  pA\ 

z  =  i~'  u  =  TT  =  Arru’  p  =  —  =  —fi2  =  —p 

Lh  U\  Ml  Pref  pc  I  M~ 

(29) 

Transformation  of  the  continuity  equations  yields  a  non- 
dimensional  form  as 


stack  level,  an  analogous  K\  must  be  determined  that  relates 
the  mass-flow  rate  through  an  individual  layer  M\  to  the 
pressure  difference  between  the  feed  manifold  and  the 
exhaust  manifold  pm  j  —  pm  e  at  that  layer. 

The  value  of  K\  follows  from  the  solution  of  the  layer 
problem.  By  definition  of  K\ 

Mi  =  Ki(p{(  0)-0)  (37) 


dwi,  ,K’M\Lh  „ 

-fiT  =  ±  pA2  iPf  -  Pc)  =  j(pf  -  pe) 


where  the  nondimensional  group  is  defined  as 


Since  it  is  assumed  that  the  exhaust  pressure  leaving  the 
(30)  layer  is  at  a  reference  value  of  zero  (See  Fig.  5),  only  the  feed 

pressure  appears  in  the  above  equation.  Upon  nondimensio- 
nalization. 


_  K’cMiLh 

‘~~pW 


(31) 


Mi 


Mi 

KipU? 


ML 

KiMi 


Ff(0) 


(38) 


Transformation  of  the  momentum  equations  leads  to  the 
following  nondimensional  form 


dgh 

dz 


Fh^hAh ,  ~~  d hh  d/3h  „ 

— -r- — Hh  =  2  uh  —  +  —  +  r  ,uh  =  0 
Mi  dz  dz 


(32) 


The  value  of  /)f(0)  is  determined  from  the  solution  of  the 
nondimensional  layer  problem.  Using  the  definition  of  Q\,  it 
follows  that 


Mh  _  Kc 

Pf{0)Mi  pf(0)Qi 


(39) 


Here,  the  nondimensional  group  is 


r. 


FhLhAh 

Mi 


u 

Mi 


(33) 


The  Hagen-Poiseuille  flow  in  the  channels  provided  a 
relationship  between  the  channel  mass-flow  rate  and  the 
pressure  drop  (Eq.  (6)).  Upon  nondimensionalization. 


mc  =  7 T  =  “  Pc)  (34) 

Mi 

The  nondimensional  mass-flow  rate  in  an  individual  channel 
(normalized  by  the  mass  flow  rate  for  the  entire  layer)  is 
proportional  to  the  nondimensional  pressure  drop  for  the 
channel.  The  nondimensional  parameter  is 


Qi 


KcMi 

Ms 


(35) 


The  relationship  between  the  continuous  and  discrete  para¬ 
meters  is 


P2\  — 


Kc 

KLh 


Q[ 


(36) 


4.  Manifold-layer  analysis 

The  foregoing  analysis  was  presented  from  the  point  of 
view  of  a  single  layer.  A  fuel-cell  stack  is  composed  of  many 
layers,  with  manifolds  supplying  and  exhausting  gas  from 
the  layers.  Fortunately,  there  is  an  analogous  analysis  for  the 
stack  of  layers.  At  the  layer  level,  the  individual  channel  flow 
is  characterized  by  the  factor  Kc  that  was  derived  from 
Hagen-Poiseuille  flow  to  relate  the  channel  mass-flow  rate 
to  the  pressure  variation  through  the  channel  (Eq.  (6)).  At  the 


Therefore,  the  value  of  K\  for  the  layer  can  be  determined 
from  the  nondimensional  solution  for  the  layer. 

With  the  value  of  Ki  in  hand,  the  analysis  for  the  stack 
follows  directly  in  the  same  form  as  the  layer  analysis.  The 
two  characteristic  nondimensional  numbers  follow  from  the 
stack  mass-flow  rate  and  the  manifold  dimensions, 


F  I  A 
AT 


2D, 


-pRef 


Fm 

AT 


(40) 


In  these  definitions,  and  throughout  the  analysis,  the  sub¬ 
script  1  for  the  layer  replaces  the  subscript  c  for  the  channel. 
Similarly,  the  subscript  m  for  the  manifold  replaces  the 
subscript  h  for  the  header.  Of  course,  there  generally  are 
a  different  number  of  layers  in  the  stack  /Vs  than  there  are 
channels  in  the  layer  /Vc. 


5.  Discretization 

Because  there  are  a  discrete  number  of  channels  that 
connect  the  feed  header  to  the  exhaust  header,  it  is  appropriate 
to  represent  and  solve  the  conservation  equations  in  discrete 
form.  Fig.  6  shows  a  header  control  volume,  which  forms  the 
structure  for  the  discrete  continuity  and  momentum  equa¬ 
tions.  For  the  purposes  of  illustration,  the  headers  are  taken  as 
rectangular,  with  height  and  width  Hi,  and  VIT .  Likewise,  the 
connecting  channels  are  rectangular,  with  dimensions  Hc  and 
Wc.  The  pitch  S  is  the  distance  between  channels  (in  the  stack 
analysis,  the  pitch  is  the  net  thickness  of  a  layer). 

There  are  Nc  discrete  channels  in  the  layer.  Consequently, 
the  headers  may  be  discretized  on  a  mesh  of  J  nodes  (Fig.  5), 
where  the  number  of  nodes  is  one  greater  than  the  number  of 
channels,  J  =  Nc+  1 .  Node  j  =  1  corresponds  to  the  left 
boundary  (z  =  0)  and  j  =  J  corresponds  to  the  right  boundary 
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Fig.  6.  Illustration  of  a  control  volume  within  a  rectangular  header, 
showing  the  entrance  to  a  rectangular  channel  connecting  the  feed  header 
to  the  exhaust  header. 


(j  =  J).  With  the  mesh  aligned  with  the  channels,  it  is  evident 
that  the  mesh  spacing  is  equal  to  the  channel  pitch,  A z  =  S. 

As  written  in  Eqs.  (8)  and  (12),  definition  of  m'  and  the 
continuity  equation  are  written  as  though  the  mass-flow  rate 
between  the  headers  varies  continuously  with  the  position  in 
the  headers  z  and  is  proportional  to  the  pressure  difference 
between  the  headers.  Accordingly,  the  proportionality  con¬ 
stant  K'q  that  relates  the  mass-flow  rate  to  the  pressure 
difference  is  written  on  a  per-unit-length  basis.  That  is, 
K'c  has  the  units  (kg-m/s-N).  However,  the  individual  chan¬ 
nels  are  spaced  along  the  headers  with  a  pitch  S  (Fig.  6). 
Therefore,  within  each  discrete  control  volume  (Fig.  5)  the 
mass-flow  rate  between  the  headers  is  determined  from  the 
Hagen-Poiseuille  flow  (Section  2.1).  The  relationship 
between  the  K'c  and  Kc  is  Kc  =  K'CS.  In  the  nondimensional 
form,  it  follows  that  the  nondimensional  group  in  the 
continuity  equation  may  be  written  as 


Q[  = 


KMh 

pAlS 


kJa_ i 
paI 


Nc  =  QiNc 


(41) 


Note  that  the  number  of  channels  and  the  nondimensional 
mesh  interval  are  related  simply  as 


Nc 


Lh  _ 

S  A  z 


(42) 


2uf 


Uf j  ~ 

As 


P f,  ~  Pfj-j 

As 


+r  \Uf  —  o. 


(2  <;</) 

(44) 


Note  that  the  convective  term  in  the  momentum  equation 
must  respect  an  upwind  difference.  Presuming  that  the 
velocity  is  positive  (i.e.  in  the  direction  of  positive  s),  the 
velocity  derivative  must  be  differenced  so  as  to  commu¬ 
nicate  information  only  in  the  flow  direction.  Failing  to 
observe  the  upwind  differencing  leads  to  numerical  instabil¬ 
ity. 

In  the  exhaust  header,  the  sense  of  the  differencing  is 
largely  opposite  that  of  the  feed  header, 

Mc;  -2  -Ply  Pc/ ) 

As  As 

=  NcQ\(ptj  -  pej),  (2  <  y  <  J)  (45) 


2m„. 

CJ 


(Mqj  Uqj_  j 

Al 


(i<y</-i) 


+  r iuej  —  o, 


(46) 


Since  the  velocity  is  also  positive  in  the  exhaust  header,  the 
upwind  differencing  of  the  convective  terms  remains  the 
same  as  in  the  feed  header. 

From  the  discrete  form  of  the  continuity  equations 
(Eqs.  (43)  and  (45))  it  is  apparent  that  the  product  NCQ\ 
can  serve  as  the  characteristic  parameter  instead  of  Qi  alone. 
Indeed,  the  solutions  presented  in  subsequent  sections  are 
characterized  in  terms  of  NCQ\,  making  the  nondimensional 
results  applicable  for  any  number  of  channels. 


5.1.  Boundary  conditions 


The  nondimensional  boundary  conditions  for  velocity  in 
the  discrete  system  are  stated  as 

Ufj  =  0,  Mc,  i  =0  (47) 


Since  the  differential  equations  are  first-order  and  because 
boundary  conditions  must  be  satisfied  at  each  end  of  the 
header  domain,  it  is  appropriate  to  forward  difference  one  of 
the  equations  and  backward  difference  the  other.  In  this  way, 
one  of  the  equations  propagates  boundary-condition  infor¬ 
mation  into  the  interior  from  each  boundary.  Even  though 
the  equations  are  coupled  and  must  be  viewed  as  a  system,  it 
can  be  helpful  to  notionally  associate  the  continuity  equation 
with  velocity  as  its  dependent  variable  and  the  momentum 
equation  with  pressure  as  its  dependent  variable.  With  this 
association  in  mind,  it  is  reasonable  to  difference  the  conti¬ 
nuity  equation  from  the  closed  header  ends  and  the  momentum 
equation  from  the  open  manifold  ends.  Specifically,  for  the 
feed  header, 

«f/+i  —  «f,  _  QiPfj  ~  Pej ) 

Az  A2 

=  -(V<A(pf.  -  Pq. ) ,  (1  <j  <  J  -  1)  (43) 


At  the  end  of  the  exhaust  header,  the  pressure  boundary 
condition  is 

Pej  =  0  (48) 

The  final  boundary  condition  is  an  implicit  one  in  which  the 
net  mass-flow  rate  for  the  layer  M\  is  specified.  The  natural 
boundary  condition  pf  \  must  be  determined  in  the  course  of 
the  solution  procedure  (discussed  in  Section  6). 

6.  Numerical  algorithm 

The  system  of  governing  equations  represents  a  bound¬ 
ary-value  problem.  Following  discretization,  the  solution  is 
accomplished  by  solving  a  system  of  nonlinear  algebraic 
equations.  The  solution  method  used  here  is  based  on  a 
modified  Newton  algorithm  [5]  and  is  implemented  in  the 
Twopnt  software  [6J. 
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The  dependent  variable  vector  0  is  comprised  of  the 
velocities  and  pressures  at  each  of  the  mesh  points  in  both 
of  the  headers.  Thus,  at  each  mesh  point  there  are  four 
dependent  variables  (1  <  k  <  4).  In  addition,  the  pressure  in 
the  feed  manifold  pQ  represents  a  single  dependent  variable. 

4*  (Pd >  “I'  •  P[\  •  Me i  •  Pe i  >  •  •  •  )  Ufj )  Pfj )  Me( ,  Fe;- 1 

•■•iMf jiPijiK’Pejf  (49) 

To  form  a  closed  system  of  algebraic  equations,  there  must 
be  a  corresponding  set  of  four  residual  equations  at  each  of 
the  mesh  points.  That  is,  there  is  a  vector  of  residuals  R  that 
is  formed  from  the  discrete  form  of  the  governing  equations 
and  boundary  conditions.  The  length  of  R  is  4/Vc  +  1. 

A  solution  vector  is  sought  for  the  solution  vector  that 
causes  the  residual  vector  to  vanish,  i.e. 

R{(j))  =  0  (50) 

The  Newton  iteration  for  solving  Eq.  (50)  may  be  written  as 

■/(0(ffl))(</>(",+1)  -  =  - R{ ( m  =  0, 1,2,3, . . .) 

(51) 


where  the  Jacobian  matrix 

J(^)  =  g  (52) 

is  evaluated  at  the  (w)th  iteration.  The  procedure  begins  with 
a  “guess”  at  the  solution  vector  After  solving  the 
system  of  linear  equations  (Eq.  (51))  for  the  “correction 
vector”  A the  (m  +  l)th  iteration  to 
the  solution  is  given  as 

^"'+0  =  +  A</>(,h)  (53) 

The  iteration  procedure  continues  until  the  correction 
becomes  negligibly  small.  Insofar  as  the  initial  iterate  lies 
within  a  domain  of  convergence,  Newton’s  method  con¬ 
verges  very  rapidly  (quadratically).  For  the  problem  at  hand, 
the  Newton  algorithm  is  quite  robust.  At  each  mesh  point, 
including  the  boundaries,  the  discrete  equations  may  be 
written  in  a  residual  form.  For  example,  the  residuals  of 
the  continuity  equation  in  the  feed  header  (A:  =  1 )  may  be 
written  as 

Rij  =  lHl"fzU'J  +NcQt(pfj  (1  <j  <  J-  1) 

RLj  =  0-  ulp  ( j  =  J ) 

(54) 

Similarly,  the  feed-header  momentum  equation  (k  =  2)  is 
written  as 


Rij  —  2ufj 


Rij  =  Po~Pfj, 


+  rmP 


(2  <j<J) 
0=  1) 

(55) 


Finally,  there  is  a  residual  equation  that  relates  the  velocity 
(from  the  feed  manifold)  into  the  feed  header  to  the  specified 
mass-flow  rate  for  the  layer.  The  residual  equation  is 

Ri  =  1  —  Mfi  (56) 

By  definition,  the  reference  velocity  U\  is  the  mean  velocity 
that  enters  the  feed  header  from  the  manifold.  Therefore,  the 
manifold  residual  equation  requires  that  u\  \  =  U\,  or 
Mf,i  =  1.  Note  that  the  associated  dependent  variable  (the 
first  variable  in  the  dependent  variable  vector)  is  the  mani¬ 
fold  pressure  p0 ,  not  a  velocity  as  might  be  expected.  This  is 
a  so-called  implicit  boundary  condition. 

The  Twopnt  boundary- value-problem  software  is 
designed  especially  to  accommodate  implicit  constraints 
or  boundary  conditions  such  as  represented  by  the  first 
residual  equation.  It  is  also  designed  to  take  advantage  of 
a  Jacobian  with  a  band  structure  such  as  in  this  problem. 


7.  Design  space 

To  explore  potential  design  space,  the  nondimensional 
problem  is  solved  for  essentially  all  plausible  combinations 
of  r i  and  NCQ\.  Fig.  7  is  a  nondimensional  map  that  provides 
quantitative  information  about  layer  performance  as  a  func¬ 
tion  of  r i  and  NCQ\.  The  flow-uniformity  contours, 

A/«c  max(mc)  —  min(mc) 

?  /  2  ,  (57) 

mc  max(mc) 

represent  the  relative  difference  between  the  mass-flow  rates 
in  each  of  the  /Vc  channels.  Assuming  that  uniform  flow  is 
desired,  a  good  layer  design  would  result  in  a  small  value  of 
A mc/mc.  The  design  space  for  Amc/mc  <0.1  is  visually 
apparent  from  Fig.  7  as  the  tent-shaped  region  below  the 


Fig.  7.  Contours  of  the  nondimensional  pressure  drop  and  channel-flow 
uniformity  as  a  function  of  and  Nc£2i. 
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10%  uniformity  contour.  Quantitatively,  this  space  is 
bounded  by 

NCQ\  <  5  x  10~2,  and  2.5TiNcQ]  <  1  (58) 

Within  the  uniform-flow  region,  the  nondimensional  pres¬ 
sure  drop  is  seen  to  depend  only  on  NcC2\  as 


A p  = 


1 


(59) 


In  the  region  above  and  to  the  right  of  the  uniformity 
contours,  some  channels  have  no  flow,  or,  potentially,  reverse 
flow.  Therefore,  because  these  regions  are  of  no  interest  to 
applications  like  fuel  cells  they  are  not  investigated. 


7.1.  Flow  distribution 


Fig.  8  illustrates  how  the  flow  distribution  in  individual 
channels  varies  for  different  combinations  of  Q\  and  Q.  The 
bar  charts,  which  are  evaluated  for  Nc  =  50  channels,  show 
the  nondimensional  mass  flow  in  each  channel.  The  inlet  is 
on  the  left  side  and  the  outlet  is  on  the  right.  By  definition  of 
the  nondimensional  mass  flow  rate  (Eq.  (34)), 

Nc 

y,  mc  =  1  (60) 

n=  1 


The  channel-flow  distribution  changes  qualitatively  in  the 
different  regions  of  nonuniform  flow.  For  low  values  of 
NCQ\,  where  the  uniformity  contours  appear  as  downward- 
sloping  lines,  the  channel-flow  distribution  is  symmetric 
with  minimum  flow  in  the  center  channels.  For  low  values 
of  r i,  where  the  uniformity  contours  appear  as  vertical  lines, 
the  channel-flow  distribution  increases  monotonically  from 
a  minimum  flow  near  the  inlet  to  a  maximum  flow  near  the 
outlet.  In  the  knee  region  where  the  uniformity  contours  turn 
the  corner  from  sloping  to  vertical  lines,  the  channel-flow 
distribution  transitions  from  the  center-low  distribution  to 
the  monotonically  increasing  distribution. 

It  is  apparent  that  the  nondimensional  pressure  drop 
increases  as  NcC2\  decreases.  At  first,  this  trend  may  seem 
counterintuitive,  since,  for  example,  lower  NCQ\  can  result 
from  lower  channel  conductance  Kc  or  lower  mass  flow  rate 
M\.  However,  note  that  the  definition  of  the  nondimensional 
pressure  can  be  written  to  involve  Q\  as 

•  2 

A/2  =  -^4  Ap  =  QiAp  (61) 

pAh  Ac 

Thus,  in  the  uniform-flow  region,  where  Ap  =  1  / NCQ\,  an 
anticipated  result  is  recovered, 


Thus,  everywhere  in  the  uniform-flow  region,  each  of  the  50  .  _  Af* 

channels  has  a  flow  rate  of  mc  =  2  x  IQ"2.  '  KCNC 


(62) 


Fig.  8.  Illustration  of  the  channel-flow  distributions  as  a  function  of  Q\  and  Fj. 


R.J.  Kee  et  at. / Journal  of  Power  Sources  109  (2002)  148-159 


157 


7.2.  Mathematical  analysis  of  the  map 


The  qualitative  behavior  of  the  solutions  can  be  under¬ 
stood  from  the  structure  of  the  differential  equations. 
Restated,  the  continuity  and  momentum  equations  are 


(63) 

=  +Q[(pf  —  pe) 

(64) 

dMf  dpf  „ 

d;  d~ 

(65) 

d Me  _  d pe 

-  i  \Mq 

d"  d" 

(66) 

Were  it  not  for  the  convective  term  in  the  momentum 
equations,  the  equations  and  boundary  conditions  preserve 
a  certain  symmetry.  For  example,  the  velocities  in  the 
channels  would  take  the  form  ue  =  1  —  iq.  The  pressure 
difference  between  feed  and  exhaust  channels  would  be 
symmetric  about  the  center  channel. 

7.2.1.  Large  T\ 

For  sufficiently  large  values  of  Fi,  the  drag  term  in  the 
momentum  equation,  Fit),  dominates  and  the  convective 
term  becomes  negligible.  Under  these  circumstances  the 
flow  must  be  symmetric.  As  illustrated  in  Fig.  8,  the  nonuni¬ 
form  channel  flow  is  always  symmetric  for  F]  >  10.  Fig.  9 
illustrates  in  more  detail  the  velocity  and  pressure  profiles  in 
the  headers  as  well  as  the  channel  mass-flow  rates  for  a 
symmetric,  but  nonuniform  flow. 

It  may  be  observed  that  Q\  and  F i  are  related  through  M\ . 
Combining  the  two  definitions  by  eliminating  M\  yields 

QiA  =  (67) 

PAh 

The  new  nondimensional  group  'F\  depends  on  the  header 
and  channel  geometry  and  the  fluid  properties,  but  not  the 
flow  rates.  Thus,  for  a  fixed  layer  and  channel  layout, 
the  product  12iFi  is  a  constant.  Taking  the  logarithm  of  the 
equation  yields 

log(Fi)  =  -log(fl1)+l0g(yi)  (68) 

which,  on  a  log-log  scale  (as  Fig.  7),  is  the  equation  of  a 
straight  line  with  a  slope  of  negative  unity.  This  observation 
is  relevant  in  understanding  the  nondimensional  mass  flow- 
rate  distribution  along  equi-uniformity  contours. 

Consider  the  governing  equations  in  the  high  F i  region, 
where  the  convective  term  is  negligible.  After  multiplying 
the  momentum  equation  by  12[,  the  system  becomes 

^-  =  ±£2[{pt-pe)  (69) 


Fig.  9.  Nondimensional  velocity  and  pressure  profiles  in  the  feed  and 
exhaust  headers  as  well  as  the  nondimensional  mass-flow  rate  in  each  of  50 
individual  channels.  At  relatively  high  r i,  the  flow  is  seen  to  be  symmetric 
but  not  uniform. 


In  this  form  it  is  clear  that  the  product  Q[p  depends  only  on 
the  nondimensional  parameter  f'J.  In  the  high  F i  region,  the 
uniformity  contours  on  the  log-log  plot  (Figs.  7  and  8)  are  a 
straight  line  with  slope  of  negative  unity.  Therefore,  along 
any  uniformity  contour,  the  product  Q\r i  =  F]  is  a  constant. 
It  follows  that  anywhere  along  a  certain  uniformity  contour 
(in  the  high  F i  region),  the  solutions  within  the  headers  must 
produce  the  same  profiles  of  Q\p(z).  Furthermore,  since  the 
nondimensional  mass-flow  rates  are 

mc  =  Q\iPi  —  Pc)  (71) 

the  nondimensional  mass-flow  rate  distribution  must  be  the 
same  anywhere  along  a  certain  uniformity  contour. 

7.2.2.  Small  T i 

For  sufficiently  small  values  of  F],  the  drag  term  in  the 
momentum  equation  becomes  negligible  and  the  convective 
term  dominates.  Consequently,  the  nonuniform  flow  must  be 
asymmetric.  Furthermore,  when  the  F \u  term  is  negligible, 
the  solutions  depend  only  on  Q\.  This  behavior  is  seen  easily 
in  Fig.  7,  where  both  the  nondimensional  pressure  and 
uniformity  contours  are  vertical  lines.  Fig.  10  illustrates 
solutions  in  the  Fi -independent  region  and  at  high  enough 
NCQ\  that  the  channel  flow  is  significantly  nonuniform. 

7.2.3.  Small  Fi  and  NCQ\ 

The  behavior  for  low  values  of  both  Q\  and  F i  can  be  seen 
by  substituting  the  continuity  equation  into  the  momentum 
equation,  yielding 


«f  [F l  -  2Q[(pf  -  pe)]  +  ^  =  0 


(72) 


Q[  ^  +  fljFifih  =  Q[  ^  +  yjfih  =  0 

d;  d- 


(70)  ue[ri  +  2Q[(pt-pe)}+^f  =  0 

d  z 


(73) 
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Fig.  10.  Nondimensional  velocity  and  pressure  profiles  in  the  feed  and 
exhaust  headers  as  well  as  the  nondimensional  mass-flow  rate  in  each  of  50 
individual  channels.  At  relatively  low  r i  and  NCQ\  near  unity,  the  flow  is 
seen  to  be  asymmetric  and  nonuniform. 


Fig.  12.  Enlarged  views  of  the  nondimensional  map  in  the  region  of  the 
knee  in  the  uniformity  contours. 


When  both  Q\  and  /)  are  sufficiently  small,  the  pressure 
gradients  are  very  small  and  the  pressure  is  nearly  uniform  in 
both  the  feed  and  exhaust  headers.  Thus,  the  channel  flow  is 
highly  uniform. 

7.2.4.  Knee  region 

In  the  region  around  T i  «  10,  where  the  drag  and  con¬ 
vective  terms  are  comparable,  the  channel  flow  transitions 
from  the  symmetric  to  the  asymmetric  patterns.  Fig.  1 1 
illustrates  a  solution  in  the  knee  region,  where  both  the 
convective  and  drag  terms  in  the  momentum  equation  are 
comparable. 

As  practical  matter,  the  best  design  space  tends  to  be  near 
the  inside  of  the  knee  of  the  uniformity  contours,  i.e.  around 


Fig.  11.  Nondimensional  velocity  and  pressure  profiles  in  the  feed  and 
exhaust  headers  as  well  as  the  nondimensional  mass-flow  rate  in  each  of  50 
individual  channels.  This  simulation  is  near  the  outside  of  the  knee  in  the 
uniformity  contours,  where  the  flow  is  in  transition  between  the  symmetric 
and  asymmetric  behavior. 


NCQ\  fts  1 0  2  and  fi  «  1.  Moving  farther  into  the  uniform 
region  usually  comes  at  the  cost  of  large  header  cross- 
sections  and  low  flow  rates,  or  large  pressure  drops  asso¬ 
ciated  with  small  channel  dimensions  (i.e.  small  Kc).  Thus, 
to  use  the  active  area  of  a  fuel-cell  layer  to  best  advantage,  a 
good  compromise  is  to  operate  near  the  regions  of  non¬ 
uniformity  with  NCQ\  as  high  as  possible.  Fig.  12,  which  is 
an  enlargement  of  the  nondimensional  map  near  the  knee  in 
uniformity  contours,  provides  more  detailed  contours  that 
facilitate  quantitative  analysis. 

8.  Channel  variations 

In  fuel-cell  applications,  the  channel  dimensions  can  be 
very  small  (order  millimeter).  Furthermore,  depending  on 
manufacturing  processes,  the  exact  dimensions  may  be 
difficult  to  control  and  they  can  have  imperfections  such 
as  in  surface  finishes.  Consequently,  it  is  interesting  to  know 
how  such  variations  may  affect  flow  characteristics. 

Consider  a  situation  in  which  the  flow  is  nominally  uni¬ 
form,  such  as  flow  represented  by  the  upper  left-hand  panel 
of  Fig.  8.  This  flow  represents  a  reasonable  design,  just 
inside  the  knee  of  the  uniformity  contours.  To  represent 
variations  in  channel  characteristics,  the  value  of  Q\  for  each 
of  50  channels  is  varied  randomly  by  ±10%  about  its 
nominal  value.  Fig.  13  illustrates  the  effect  on  channel  flow 
distribution. 

From  Eq.  (7),  it  is  seen  that  the  channel  conductance  Kc 
depends  strongly  on  the  channel-cross-section  dimensions. 
Specifically,  Kc  depends  on  the  cube  of  a  characteristic 
channel  dimension  (e.g.  its  height).  Therefore,  a  10%  varia¬ 
tion  in  channel  height  causes  a  much  larger  variation  in  Kc, 
and  consequently  Q\.  As  a  result,  to  assure  uniform  flow. 
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Fig.  13.  Nondimensional  velocity  and  pressure  profiles  in  the  feed  and 
exhaust  headers  as  well  as  the  nondimensional  mass-flow  rate  in  each  of  50 
individual  channels.  This  simulation  is  near  the  inside  of  the  knee  in  the 
uniformity  contours,  where  the  flow  is  nearly  uniform.  To  represent  the 
effect  of  fabrication  tolerances,  the  value  of  Q\  is  varied  randomly  by 
±10%  in  each  channel. 

there  is  a  very  important  demand  on  fabrication  processes  to 
maintain  control  of  channel  dimensions. 


9.  Conclusions 

An  analysis  has  been  developed  to  characterize  the  flow 
distribution  in  planar  fuel-cell  layers  and  stacks.  Using 
nondimensional  variables,  the  results  are  presented  in  gra¬ 
phical  form  that  can  be  applied  generally.  While  the  analysis 


is  principally  motivated  by  fuel-cell  applications,  other 
applications  like  plate-and-frame  heat  exchangers  can  ben¬ 
efit  from  the  results. 
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